Entanglement-assisted quantum low-density parity-check codes 
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This paper develops a general method for constructing entanglement-assisted quantum low-density 
parity-check (LDPC) codes, which is based on combinatorial design theory. Explicit constructions 
are given for entanglement-assisted quantum error-correcting codes (EAQECCs) with many desir- 
able properties. These properties include the requirement of only one initial entanglement bit, high 
error correction performance, high rates, and low decoding complexity. The proposed method pro- 
duces several infinite families of new codes with a wide variety of parameters and entanglement 
requirements. Our framework encompasses the previously known entanglement-assisted quantum 
LDPC codes having the best error correction performance and many new codes with better block 
error rates in simulations over the depolarizing channel. We also determine important parameters 
of several well-known classes of quantum and classical LDPC codes for previously unsettled cases. 
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I. INTRODUCTION 

This paper develops a general combinatorial method 
for constructing quantum low-density parity-check 
(LDPC) codes under the entanglement-assisted stabi- 
lizer formalism established by Brun, Devetak, and Hsich 
Our results include many new explicit construc- 
tions for entanglement-assisted quantum error-correcting 
codes for a wide range of parameters. We also prove a 
variety of new results for classical error-correcting codes, 
which directly apply to the quantum setting. Most of 
the quantum codes designed in this paper achieve high 
error correction performance and high rates while requir- 
ing prescribed amounts of entanglement. These codes 
can be efficiently decoded by message-passing algorithms 
such as the sum-product algorithm (for details of itera- 
tive probabilistic decoding, see [2(). 

The existence of quantum error-correcting codes was 
one of the most important discoveries in quantum infor- 
mation science [H, Q • Unfortunately, most of the known 
quantum error-correcting codes lack practical decoding 
algorithms. 

In this paper, we focus on the use of LDPC codes in 
a quantum setting. Classical LDPC codes @ can be ef- 
ficiently decoded while achieving information rates close 
to the classical Shannon limit This extends to the 

quantum setting: the pioneering works of Hagiwara and 
Imai [9| and MacKay, Mitchison, and McFadden [lOj pre- 
sented quantum LDPC codes which surpassed, in simu- 
lations, all previously known quantum error-correcting 
codes. Their quantum codes have nearly as low decoding 
complexity as their classical counterparts. 

However, most of the previous results concerning quan- 
tum LDPC codes and related efficiently decodable codes 
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have relied on the stabilizer formalism, which severely 
restricts the classical codes which can be used. The diffi- 
culty in developing constructions for non-stabilizer codes 
was also a substantial obstacle. 

Our results use the newly developed theory of 
entanglement-assisted quantum error-correcting codes 
(EAQECCs) HUMl. The entanglement-assisted stabi- 
lizer formalism allows the use of arbitrary classical binary 
or quaternary linear codes for quantum data transmis- 
sion and error correction by using shared entanglement 
[Til [l5| . Previous work related to entanglement-assisted 
quantum LDPC codes is due to Hsieh, Brun, and Deve- 
tak PJ] and Hsieh, Yen, and Hsu [l7|. 

The major difficulty in using classical LDPC codes in 
the entanglement-assisted quantum setting is that very 
little is known about methods for designing EAQECCs 
requiring desirable amounts of entanglement. While 
entanglement-assisted quantum LDPC codes can achieve 
both notable error correction performance and low de- 
coding complexity, the resulting quantum codes might 
require too much entanglement to be usable; in general 
entanglement is a valuable resource 15]. In some situ- 
ations, one might wish to effectively take advantage of 
high performance codes requiring a larger amount of en- 
tanglement p],[l2|. To the best of the authors' knowledge, 
no general methods have been developed which allow the 
code designer flexibility in choice of parameters and re- 
quired amounts of entanglement. 

Our primary focus in this paper is to show that it is 
possible to create infinite classes of EAQECCs which con- 
sume prescribed amounts of entanglement and achieve 
good error correction performance while allowing efficient 
decoding. Our methods are flexible and address various 
situations, including the extreme case when an EAQECC 
requires only one preexisting entanglement bit. 

The entanglement-assisted quantum LDPC codes 
which we construct include quantum analogues of the 
well-known finite geometry LDPC codes originally pro- 
posed by Kou, Lin, and Fossorier [lj| (see also [la. |20|). 



2 



and LDPC codes from balanced incomplete block de- 
signs that achieve the upper bound on the rate for a 
classical regular LDPC code with girth six proposed in- 
dependently by several authors (see [2lJ and references 
therein). Some classes of our codes outperform previ- 
ously proposed quantum LDPC codes ha ving the best 
known error correction performance 0, [ToL Il6l Il7| . 

Our primary tools come from combinatorial design the- 
ory, which plays an important role in classical coding the- 
ory [22j and also gave several classes of stabilizer codes 
in quantum coding theory (23l - [27| . The use of combi- 
natorial design theory allows us to exactly determine or 
give tighter bounds on the parameters of the finite geom- 
etry LDPC codes in both quantum and classical settings. 
Comprehensive lists of the parameters of these codes are 
given in Tables IXlVl and IXVl in Appendix |B] 

In Section lU we outline our framework for design- 
ing entanglement-assisted quantum LDPC codes by us- 
ing combinatorial design theory. Section IIIII gives ex- 
plicit constructions for entanglement-assisted quantum 
LDPC codes based on finite geometries and related com- 
binatorial structures. New results concerning the well- 
known classical finite geometry LDPC codes are also 
given in this section. Section IIVI presents simulation re- 
sults of our entanglement-assisted quantum LDPC codes 
and discusses their performance over the depolarizing 
channel. Section [V] contains concluding remarks and dis- 
cusses some related problems that can be treated with 
the techniques developed in this paper. 



II. COMBINATORIAL 
ENTANGLEMENT-ASSISTED QUANTUM LDPC 
CODES 

In this section we give a general construction method 
for entanglement-assisted quantum LDPC codes based 
on combinatorial designs. We do not describe the theory 
of classical LDPC codes in detail here, instead referring 
the reader to 0, HH and references therein. Relations be- 
tween quantum error-correcting codes and LDPC codes 
are concisely yet thoroughly explained in [l(| fl6| - Ba- 
sic notions related to LDPC codes and their relations to 
combinatorial designs can be found in [28j . For a de- 
tailed treatment of the entanglement-assisted stabilizer 
formalism, we refer the reader to 0, [l2l - [l4| . 

In Subsection lll Al we introduce necessary notions from 
coding theory and combinatorial design theory. A gen- 
eral method for designing entanglement- assisted quan- 
tum LDPC codes is presented in Subsection lll Bl 



A. Preliminaries 

An [[n, fc; c\] entanglement- assisted quantum error- 
correcting code (EAQECC) encodes k logical qubits into 
n physical qubits with the help of c copies of maxi- 
mally entangled states. As in classical coding theory, 



n is the length of the EAQECC, and k the dimen- 
sion. We say that the EAQECC requires c ebits. An 
[[n, k; c]] EAQECC with distance d will be referred to as 
an [[n, fc, d; c]] code. 

The rate of an [[n,fc;c]] EAQECC is defined to be |. 
The ratio is called the net rate. The latter figure de- 
scribes the rate of an EAQECC when used as a catalytic 
quantum error-correcting codes to create c new bits of 
shared entanglement [l|, Il2| . 

Throughout this paper, matrix operations are per- 
formed over F2, the finite field of order two. The ranks 
of matrices are also calculated over F 2 . 

We employ the Calderbank-Shor-Steane (CSS) con- 
struction 0, H G3, HI]. Usually the CSS construction 
uses a minimal set of independent generators to construct 
an EAQECC. Hence, the construction is often described 
by using a classical binary linear code with a parity- 
check matrix of full rank. However, in actual decoding 
steps, sparse-graph codes may take advantage of redun- 
dant parity-check equations to improve error correction 
performance. Because the extended syndrome can be ob- 
tained in polynomial time without additional quantum 
interactions, we use the following formulation of the CSS 
construction for EAQECCs. 

Theorem 1 (Hsieh, Brun, and Devetak [16]) // 

there exists a classical binary [n, k,d] code with parity- 
check matrix H , then there exists an [[n, 2k — n + c, d; c]] 
EAQECC, where c = rank HH T . 

Note that H may contain redundant rows which are re- 
lated only to classical operations to infer the noise by a 
message-passing algorithm. 

We apply Theorem Q] to classical sparse-graph codes. 
An LDPC code is typically defined as a binary linear 
code with parity-check matrix H in which every row and 
column is sparse. In this paper we consider LDPC codes 
with parity-check matrices whose rows and columns con- 
tain only small numbers of ones so that simple message- 
passing algorithms can efficiently give good performance 
in decoding. 

Proposition 2 An LDPC code with parity-check matrix 
H with n columns and minimum distance d defines a 
classical binary [n,n— rankTJ, d] code, which yields an 
[[n,n-2rankiJ+ mnkHH T ,d] r&nkHH T }} EAQECC. 

The Tanner graph of an m x n parity-check matrix H 
is the bipartite graph consisting of n bit vertices and m 
parity-check vertices, where an edge joins a bit vertex to 
a parity-check vertex if that bit is included in the cor- 
responding parity-check equation. A cycle in a graph is 
a sequence of connected vertices which starts and ends 
at the same vertex in the graph and contains no other 
vertices more than once. The girth of a parity-check ma- 
trix is the length of a shortest cycle in the corresponding 
Tanner graph. Short cycles can severely reduce the per- 
formance of an otherwise well-designed LDPC code. In 
fact, one of the greatest obstacles to the development of 
a general theory of LDPC codes in the quantum setting 
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is the difficulty of avoiding cycles of length four (See, for 
example, 0, [T3, [3(3, HU). In order to improve error cor- 
rection performance, we generally only treat LDPC codes 
with girth at least six. 

The weight of a row or column of a binary matrix 
is its Hamming weight, that is, the number of ones in 
it. An LDPC code is regular if its parity-check matrix 
H has constant row and column weights, and irregular 
otherwise. Regular LDPC codes are known to be able 
to achieve high error correction performance. Irregular 
LDPC codes allow the code designer to optimize charac- 
teristics of performance by a careful choice of row weights 
and column weights 0-11]. 

We now define several combinatorial structures, which 
we will need in Subsection III Bl and the subsequent sec- 
tions. For additional facts and desi gn t heoretical results, 
the interested reader is referred to |32j |. 

An incidence structure is an ordered pair (V, B) such 
that V is a finite set of points, and B is a family of subsets 
of V , called blocks. A point-by-block incidence matrix of 
an incidence structure (V, B) is a binary v x b matrix 
H = (h-ij) in which rows are indexed by points, columns 
are indexed by blocks, and hij — 1 if the ith point is 
contained in the jth block, and hij = otherwise. A 
block-by-point incidence matrix of (V, B) is the transposed 
point-by-block incidence matrix H T . 

Any LDPC code can be associated with an incidence 
structure by interpreting its parity-check matrix as an 
incidence matrix. The converse also holds as long as the 
considered incidence matrix is sparse. 

The current paper will focus on incidence struc- 
tures which have been extensively studied in combina- 
torics. This allows us to effectively exploit combinato- 
rial design theory to develop a framework for designing 
entanglement-assisted quantum LDPC codes. 

A 2-(v,fj,,X) design is an incidence structure (V,B), 
where V is a set of cardinality v and B is a family of 
/i-subsets of V such that each pair of points is contained 
in exactly A blocks. We will refer to the parameters v, ll, 
and A as the order, block size, and index of a 2-design. 
Note that the block size of a 2-design is usually written as 
k in the combinatorial literature. To avoid any confusion 
with the dimension of a code, we use ll instead. 

The number b = \B\ of blocks in a 2-(v, fx, A) design is 
determined by the design parameters: 



(1) 



fl(fj, - 1) 

A 2-design is called symmetric if b — v. 

Every point of a 2-{v,ll,X) design occurs in exactly r 
blocks, where 



r = -A. 



Ll-1 



(2) 



The number r is called the replication number of the de- 
sign. A point-by-block incidence matrix H of a 2-(v, ll, A) 
design satisfies the equation 



where / is the identity matrix and J is the v x v all-one 
matrix. Since r and b are integers, it follows that the 
following two conditions 



A(v-1) = (mod/i-1), 
Xv(v - 1) = (mod fi(n - 1)) 



(4) 



are necessary conditions for the existence of a 2-(v, ll, A) 
design. 

If the block size li and index A are relatively small, an 
incidence matrix of a 2-(v, ll, A) design is sparse. Hence, a 
point-by-block incidence matrix of a 2-{v, ll, A) design can 
be viewed as a parity-check matrix H of a regular LDPC 
code with constant row weight r and constant column 
weight ll. Similarly, a block-by-point incidence matrix 
defines a code with constant row weight ll and constant 
column weight r. In this paper, incidence matrices will 
generally be point-by-block unless it is specifically noted 
otherwise. In the cases when block-by-point matrices are 
desirable, the notation H T will be used. 

A substantial part of this paper deals with one of the 
most fundamental incidence structures in combinatorial 
design theory. A Steiner 2-design, denoted by S{2,ll,v), 
is a 2-(v, ll, 1) design. A Steiner triple system of order 
v, denoted by STS(w), is a Steiner 2-design with block 
size three. The S(2, ll, v)s are trivial Steiner 2-designs if 
V < /i. We generally do not consider trivial designs to be 
Steiner 2-designs unless they play an important role. 

It is easy to see that both point-by-block and block- 
by-point incidence matrices of an S(2,ll, v) give regular 
LDPC codes with girth six (see, for example, [33j). 



B. General combinatorial constructions 

In this subsection we present a general framework for 
designing entanglement-assisted quantum LDPC codes 
based on combinatorial design theory. Specialized con- 
struction methods for desirable EAQECCs in this frame- 
work will be illustrated in Section IHIl 

The following propositions are derived from Theorem 
[T]by using incidence matrices as parity-check matrices of 
binary LDPC codes. 

Proposition 3 Let H be a point-by-block incidence ma- 
trix of an incidence structure (V,B). Then there exists a 
[[\B\,\B\-2T&nkH+ rank HH T ; rankHH T }} EAQECC. 

Proposition 4 Let H T be a block-by-point incidence 
matrix of an incidence structure (V,B). Then there 
exists a [[\V\, \V\ - 2rank# + rankff T i/; rankff T iJ]] 
EAQECC. 

We employ the following two theorems. 

Theorem 5 (Hillebrandt [34]) The rank of an inci- 
dence matrix H of an S(2,ll,v) satisfies the following 
inequalities: 



HH T = (r - X)L + A J, 



(3) 




1 (v-l^y-v) 
4 H 



< ranki/ < v. 
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Theorem 6 (Hamada |35l |) If H is an incidence ma- 
trix of an S(2, /x, v) with even replication number r = 
then 



rank H = 



v — 1 when /i is even, 
v or v — 1 when [i is odd. 



We now give three simple constructions by applying 
Propositions [3] and |4] to incidence matrices of Steiner 
2-designs. These constructions will be specialized and 
modified to give desirable codes. 

Theorem 7 (High-Rate 1-Ebit Code) Let H be a 

point-by-block incidence matrix of an S(2,/j,,v). Suppose 
r = * s odd. Then H has row weight r, column weight 
H, girth 6, and the corresponding [[n, k; c]] EAQECC sat- 
isfies the following conditions: 



v{v — I) 
Ufa - I)' 



vr „ vr „ 
2v+l <k< 2 




1 (v-l^y-u) 
4 H 



+ 1, 



c = f. 

Proof. By Proposition [3] and Theorem [5l it suffices to 
prove that rank HH T = 1. Because r is odd, Equation 
(13]) reduces to HH T = J, which implies that the rank of 
HH T is equal to one. □ 

Theorem 8 (High-Rate High-Consumption Code) 

Let H be a point-by-block incidence matrix of an 
S(2,/i,v). Suppose r = ^Ej * s even. Then H has row 
weight r, column weight fi, girth 6, and the corresponding 
[[n, k; c]] EAQECC satisfies the following conditions: 



v(v — 1) 
fi(n - 1)' 



k = 



I' 



+ f when fi is even, 

1 or — — v — 1 when a is odd, 



Proof. By Proposition [3] and Theorem [6l it suffices to 
prove that rank HH T = v — 1 . Because r is even, Equa- 
tion ((3]) reduces to 



HH 1 = 



f 
f 



f f ••• 



that is, a matrix containing zeros on the diagonal and 
ones in the other entries. Because r — ' 1 

odd. Hence, we have rank HH T = v 



is even, v is 
□ 



1 as desired. 



Theorem 9 (Low-Rate High-Redundancy Code) 

Let H T be a block-by-point incidence matrix of an 
S(2,/j,,v). Then H has row weight fi, column weight 
r, girth 6, and the corresponding [[n, fc; c\] EAQECC 
satisfies the following conditions: 



k < v - 2 




1 (v-l)(v-H) 
4 



c, 



c > f. 



Proof. Let H T be a block-by-point incidence matrix 
of an S(2,fi,v). Since any non-trivial S(2,/i,v) contains 
a pair of blocks that share exactly one point, we have 
rankiJ T iJ > f. Applying Proposition [4] to Theorem [5] 
completes the proof. □ 

It is worth mentioning that a weaker version of The- 
orem [7] was used in the context of integrated optics and 
photonic crystal technology [25]. Also notable is that 
Theorems[7]and[5]can be easily extended to the case when 
preexisting entanglement is not available. For example, 
quantum LDPC codes that do not require entanglement 
can be obtained by applyin g th e extra column method 
used in Construction U in [l(| and the CSS construc- 
tion to S(2,[i,v)s in the same manner as in Proposition 
[3] Aly's construction for quantum LDPC codes [23| is a 
special case of this extended method. Djordjevic's con- 
struction for quantum LDPC codes 24| can be obtained 
by applying the CSS construction to 2-designs of even 
index in the same way as in Proposition [3] 

The existence of 2-designs is discussed in Appendix [Al 
which provides Steiner 2-designs necessary to obtain sev- 
eral infinite families of new entanglement-assisted quan- 
tum LDPC codes from Theorems [3 [H and [9l Before 
applying our theorems to specific S(2,fi,v)s, we explore 
general characteristics of our EAQECCs and further de- 
velop methods for designing desirable codes. 

Theorem [7] yields entanglement-assisted quantum 
LDPC codes with very high net rates and various lengths 
while requiring only one ebit. Theorem [S] gives codes 
which have very high net rates and naturally take ad- 
vantage of larger numbers of cbits when there is an 
adequate supply of entanglement. Since rankHH T < 
rankH holds for any parity-check matrix H, the required 
amounts of entanglement of high rate codes in Theorem 
[5] are expected to be relatively low when compared with 
randomly chosen codes of the same lengths. Theorem 
[5] generates entanglement-assisted quantum LDPC codes 
which can correct many quantum errors by taking advan- 
tage of the higher redundancy. The high error correction 
performance of these codes will be demonstrated in sim- 
ulations in Section HVl 

When a parity-check matrix H of an S(2,fj,,v) is of 
full rank v, the corresponding classical LDPC code in 
Theorems [7J and [5] achieves an upper bound on the rate 
for an LDPC code with girth six. 
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Theorem 10 (MacKay and Davey (36[) Let H be a 

v x n parity-check matrix of a classical regular LDPC 
code of length n, column weight fi, and girth 6. Let also 
ranki/ = v. Then it holds that n < v ^ v ~ 1 ,\ , where equal- 
ity holds if and only if H is an incidence matrix of an 
5(2,,i,«). 

It follows that EAQECCs based on Steiner 2-designs 
achieve the highest possible net rates for quantum LDPC 
codes with girth at least six constructed from full 
rank parity-check matrices with constant column weights 
through the CSS construction. 

The rank of an incidence matrix of an S(2,fi,v) may 
not be full depending on the structure of the design. 
If one wishes a parity-check matrix to be regular and 
full rank at the same time, it is important to choose 
an 5(2, /j,, v) with a full rank incidence matrix. This 
can always be done for the case when fi = 3 except for 
v = 7 [37[ • For a more detailed treatment of the ranks of 
5(2, /z, v)s, we refer the reader to [H, 155, l39j. 

In general, the code minimum distance plays less of 
a role in the performance of sum-product decoding than 
maximum likelihood decoding [Toj . Therefore, we explore 
in detail the distance d of [[n, k, d; c]] EAQECCs based on 
LDPC codes only when it is of great theoretical interest. 
Because codes derived from finite geometries are of great 
importance in coding theory, the distances of EAQECCs 
obtained from finite geometries will be investigated in 
detail in Section IIIII 

Here we briefly review the minimum distances of LDPC 
codes based on Steiner 2-designs. A pair of 5(2, /i,«)s 
which are not mutually isomorphic may give different 
minimum distances. The tightest known upper and lower 
bounds on the minimum distance of an LDPC code based 
on an STS(u) can be found in the very large scale inte- 
gration (VLSI) literature as bounds on even-freeness. 



Theorem 11 (Fujiwara and Colbourn [40]) The 

minimum distance d of a classical binary linear code 
whose parity- check matrix forms an incidence matrix of 
a non-trivial STS(v) satisfies 4 < d < 8. 

A carefully chosen triple system can have a good topo- 
logical structure which gives good decoding performance. 
If conditions require larger minimum distances, the code 
designer may use either block-by-point incidence matri- 
ces, or 5(2, /i, v)s of larger block sizes. For known results 
on minimum distances, girths, and related characteris- 
tics of LDPC codes based on combinatorial designs, the 
reader is referred to (40l - |42l | and references therein. 

In what follows, we describe general guidelines for 
designing entanglement-assisted quantum LDPC codes 
with desired parameters and properties by exploiting 
codes we have presented in this section. 

We first consider an [[n, k; c]] EAQECC requiring only 
a small amount of entanglement. The extreme case is 
when c = 1. The following theorem gives infinitely many 
such EAQECCs having extremely high rates and low de- 
coding complexity. 



Theorem 12 Let v and /j, be positive integers satisfying 
v — 1 = (mod fi — 1) and v(v — 1) = (mod /i(/i — 1)). 
Suppose also that * s odd. Then for all sufficiently 
large v and some k satisfying the condition of Theorem 



[7] there exists an 



[$E$,A;1]] EAQECC. 



Proof. Use Theorem |4"TI from Appendix A, which guar- 
antees the existence of an 5(2, fj,,v) for all sufficiently 
large v, and apply Theorem [7] □ 
Similarly, applying Theorem [7] to known 5(2, /i, v)s 
with small v discussed in Appendix [X] gives [[n, k; 1]] 
EAQECCs of shorter length n. 

In general, the error floor of a well-designed LDPC 
code is not dominated by low-weight codewords. 
Nonetheless, it is desirable to carefully choose an 
5(2, /z, v) when applying our simple constructions so that 
the resulting code has a promising topological structure. 
While incidence matrices of 5(2, /i, v)s have long been in- 
vestigated in various fields, it appears to be difficult to 
achieve the known upper bounds on the minimum dis- 
tance of an LDPC code based on an incidence matrix of 
an 5(2, v). In fact, it is conjectured that the known 
upper bounds are generally not achievable even for the 
3 ~ 



41] 



case /j, ■■ 

An STS is 4- even- free (or anti-Pasch) if its incidence 
matrix gives a classical LDPC code with minimum dis- 
tance five or greater. A 4-even-free STS(i;) exists for all 
v ^ 7, 13 satisfying the necessary conditions (g]) [43T ]. It 
is conjectured that an incidence matrix of a 4-even-free 
STS(w) gives the largest possible minimum distance [ilj ]. 

Theorem 13 There exists a [[^~^-,k,d; 1]] EAQECC 

with k > 1 - ) — 2v + I and d > 5 for every v = 3, 7 
(mod 12) except for v = 7. 

Proof. If v = 3, 7 (mod 12), then the replication num- 
ber of an STS(w) is odd. Applying Theorem [7] to a 4- 
even-free STS(w) completes the proof. □ 
A block-by-point incidence matrix of a symmetric 
5(2, /i,w) can also be viewed as a point-by-block inci- 
dence matrix of a Steiner 2-design of the same parameters 
LIT] . Hence, Theorems [7] and [9] can overlap when sym- 
metric designs are employed. This special case gives the 
EAQECCs with c — 1 and good error correction perfor- 
mance originally presented in [l7^ . For completeness, we 
give a simple proof by using the following two theorems. 



Theorem 14 For every integer t > 1 there exists a sym- 
metric 5(2, 2* + 1, 4* + 2* + 1) whose incidence matrix H 
satisfies rank H = 3 + 1. 

Proof. Take as 5(2, 2* + 1, 4* + 2* + 1) the Desarguesian 
projective plane of order 2*, whose incidence matrix has 
rank 3* + 10. □ 



Theorem 15 (Calkin, Key, and de Resmini [45]) 

Let H T be a block-by-point incidence matrix of a sym- 
metric 5(2,2* + 1 , 4* + 2* + 1) being the Desarguesian 
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projective plane PG(2,2 t ). Then H T defines a classical 
binary linear [4* + 2* + 1, 4* + 2* - 3*, 2* + 2] code. 

Now as a corollary of Theorems [7J and O and the pre- 
ceding two theorems, we obtain the following result. 

Theorem 16 For every integer t > 1 there exists a [[4* + 
2* + 1, 4* + 2* - 2 • 3*, 2* + 2; 1]] EAQECC. 

EAQECCs of this kind can be seen as quantum ana- 
logues of special Type I PG-LDPC codes, which have 
notable error correction performance in the classical set- 
ting [l8l - [2"(il ] . Because of the direct correspondence be- 
tween entanglement-assisted quantum codes and classical 
codes, these EAQECCs inherit excellent error correction 
performance while consuming only one initial ebit. We 
will further investigate entanglement-assisted quantum 
LDPC codes based on S(2,n,v)s with large minimum 
distances in Section ITO1 

Next we present general combinatorial methods for de- 
signing EAQECCs with relatively small c and better er- 
ror correction performance. The main idea is that we 
discard some columns from an incidence matrix of an 
5(2, (j., v) and then apply Proposition^ as we did in The- 
orem [Jj Our methods encompass the rate control tech- 
nique for classical LDPC codes proposed in [46| as a spe- 
cial case. 

Let (V,B) be an 5(2, fi,v). Take two subsets V C V 
and B 1 C B. The pair (V',B') is called a proper subde- 
sign of block size \i if it is an 5(2, \i,\V'\). Since we do 
not consider other kinds of subdesigns, we simply call a 
proper subdesign (V',B') of block size p, a subdesign. A 
pair of subdesigns (V',B') and {V",B") of an S(2,fj,,v) 
are point-wise disjoint if V PI V" = 0. 

Theorem 17 Let (V,B) be an S(2,fi,v) with odd r = 
^Ej- Assume that (V,B) contains j point-wise mutu- 
ally disjoint subdesigns (Vi,Bi), 1 < i < j, such that 
[J J i= i Vi C V and each (Vi,Bi) has odd replication num- 
ber. Then there exists an [[n,k;c]] EAQECC satisfying 
the following conditions: 



v(v — 1) 



c = j + l. 

Proof. Take an arbitrary incidence matrix H of an 
5(2, (j,,v) with odd r. Delete j point-wise mutually dis- 
joint subdesigns (V i: Bi) each of which has odd replica- 
tion number. It is always possible to reorder the rows 
and columns of the resulting incidence matrix H' such 
that H'H' T has the form: 



H'H' T = 



J J 
J Oi 

J J 



J 
J 



where 0^ is a \Vi\ x \Vi\ zero matrix and each J is an 
all-one matrix of appropriate size. It is easy to see that 
rank H'H lT = j + 1. Applying Proposition [3] to H' com- 
pletes the proof. □ 

Deleting subdesigns always shortens the length of the 
corresponding code. Discarding columns will not de- 
crease the minimum distance or the girth. The rank of 
the parity-check matrix is unlikely to change. In this 
sense, we expect EAQECCs obtained through subdesign 
deletion to have better error correction performance than 
the original code. We will demonstrate this effect in sim- 
ulations in Section HVl 

In general, deleting a subdesign makes a parity-check 
matrix slightly irregular. If this irregularity is not desir- 
able because of particular circumstances or conditions, it 
can be alleviated by discarding more point-wise disjoint 
subdesigns. In fact, if we delete subdesigns of the same 
order such that each point belongs to one deleted subde- 
sign, we obtain a regular parity-check matrix again. The 
following construction demonstrates this. 

Let (V, B) be an S(2,fi,v) and S a set of Steiner 2- 
designs (Vi,Bi), 1 < i < \S\, where V\, . . . , Vjs| partition 
V, that is, (j Vi = V and V l r\V j = for all i ^ j. Then S 
is called a Steiner spread in (V, B) if each (Vi,Bi) forms 
a subdesign 5(2, fi, \V t \) of (V,B). 

Theorem 18 Let (V,B) be an 5(2, fx,v) with odd repli- 
cation number r — ^Ej- Assume that (V,B) contains 
a Steiner spread S, where each subdesign (Vi,B{) has 
odd replication number. Then there exists an [[n, fc; c]] 
EAQECC satisfying the following conditions: 



v(v — 1) 
MO - !) 



|5| — 1 when \S\ is odd, 
\S\ when 151 is even. 



Moreover, if \Vi\ — \Vi'\ — w for all i and i' , then the 
parity-check matrix of the corresponding LDPC code is 
regular and has row weight r — ^-f and column weight 
fi. 

Proof. Let H be an incidence matrix of an 5(2, fi,v) 
with odd r which contains a Steiner spread S. Delete all 
members of the Steiner spread from (V, B). By following 
the same argument as in the proof of Theorem 1171 it is 
straightforward to see that rank HH T = \S\ — 1 when 
|5| is odd, and |<S| otherwise. If \Vi\ = \Vi'\ = w for all i 
and i' , each subdesign has the same replication number 
^Ej- Hence, the resulting code is regular. □ 
When there is an adequate supply of entanglement, it 
may be acceptable to exploit a relatively large amount 
of entanglement to improve error correction performance 
while keeping similar characteristics of high rate codes. 
Deleting an 5(2, /i, w) with even replication number 
increases the required amount of entanglement to a 
slightly larger extent. 
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Theorem 19 Let (V,B) be an S(2,/j,,v) with odd repli- 
cation number r = -^y- Assume that (V,B) contains j 
point-wise mutually disjoint subdesigns (Vi,Bi), 1 < i < 
j, such that IB =1 Vi C V and each (V^, £?,-) /ias even repZ«- 
cation number. Then there exists an [[n, k; c]] EAQECC 
satisfying the following conditions: 



v(v — 1) 
fi(H - 1) 



i=l 



+ 



Moreover, if the subdesigns (Vi, Bi) for 1 < i < j /orm a 
Steiner spread with \Vi\ = \Vi'\ = w for all i and i' , then 
the parity-check matrix of the corresponding LDPC code 
is regular and has row weight r — an d column weight 
fjt. 

Proof. Take an arbitrary incidence matrix H of an 
S(2,fi,v) with odd r. Delete j point-wise mutually dis- 
joint subdesigns (Vi,Bi) each of which has even replica- 
tion number. If U<=i Vi £ V, ^ is always possible to 
reorder the columns of the resulting incidence matrix H' 
such that H'H lT is of the form: 



H'H' T = 



J J 
J h 

J J 



J 
J 



where h is the \Vi\ x \Vi\ identity matrix and each J is an 
all-one matrix of appropriate size. Because each Ii has 
Vi independent rows and each \ Vi\ is odd, jankH' H lT = 
J2i=i (\ v i \ - 1) + 1< Applying Proposition [3] to H' gives 

c = ELi (\ V i\ ~ !) + L If ULi Vi = V, we have identity 
matrices across the diagonal of H 1 H' T . Hence, we have 
c = J2l=i (\Vi \ — 1) + 1 again. If each Vi is of the same 
size, it is straightforward to see that the resulting code 
is regular. □ 

When irregularity in a parity-check matrix is accept- 
able or favorable, the code designer can combine the tech- 
niques of Theorems [TTl UHl and [HI The required amount 
of entanglement is readily computed by the same argu- 
ment as above. 

In general, subdesign deletion changes the parameters 
of a code in a gradual manner. Hence, these techniques 
are also useful when one would like an EAQECC of spe- 
cific length or dimension. While we only employed The- 
orem [7] in the above arguments, Theorem [5] can also be 
used in a straightforward manner to fine-tune the param- 
eters of EAQECCs. 

In order to exploit the subdesign deletion techniques, 
one needs Steiner 2-designs having subdesigns or prefer- 
ably Steiner spreads of appropriate sizes. We conclude 
this section with a brief review of known general results 



and useful theorems for finding 5(2, fi, v) with subdesigns 
and Steiner spreads. For a more thorough treatment, the 
reader is referred to [lH, \^ and references therein. 

The well-known Doyen- Wilson theorem [48[ states that 
one can always find an STS(u) containing an STS(w;) as 
a subdesign as long as both v and w satisfy the necessary 
conditions for the existence of an STS and v > 2w + 1. 
The following is a general asymptotic theorem on Steiner 
2-designs having subdesigns. 

Theorem 20 (Fujiwara (49fl ) Let /i > 2 be a positive 
integer and w = 1 (mod fj,(fx — 1)). Then there exist 
a constant number wq depending on fi, and a constant 
number vq depending on w and fi such that if w > wq 
and v > vq satisfies the conditions v — 1 = (mod fj, — 1) 
and v(v — 1) = (mod fi(fj, — 1)), then there exists an 
5(2, n,v) having an S(2,n,w) as a subdesign. 

Theorem [20l states that one can always find an 5(2, fi, v) 
having an 5(2, fi, w) as a subdesign as long as v is a suf- 
ficiently large integer satisfying the necessary conditions 
(U|) and w is a sufficiently large integer satisfying w = 1 
(mod fi(fi — 1)). 

Steiner spreads are closely related to a special kind of 
combinatorial design. A group divisible design (GDD) 
with index one is a triple (V, Q , B), where 

(i) V is a finite set of elements called points, 

(ii) Q is a family of subsets of V , called groups, which 
partition V, 

(iii) B is a collection of subsets of V, called blocks, such 
that every pair of points from distinct groups occurs 
in exactly one block, 

(iv) \G n B\ < 1 for all G £ Q and B e B. 

If all groups are of the same size g, all blocks are of the 
same size fi, and \Q\ = t, one refers to the design as a 
^-GDD of type gK 

Theorem 21 The existence of an S(2,fi,g) and a fi- 
GDD ( V, Q, B) of type g f with index one implies the exis- 
tence of an S(2,fi,gt) having a Steiner spread S, where 
each member of S is an S(2,fi,g). 

Proof. Let (V, Q, B) be a ^-GDD of type g l with in- 
dex one and (V ,B') an S(2,fi,g). For each G € Q, we 
construct an S(2,fi,g), (G,B' G ), by mapping each point 
of (V, B') to an element of G by an arbitrary bijection 
ttg : V — > G. Define C = UggG &g- ^ ^ s straightforward 
to check that (V, BUC) is an 5(2, fi, gt) having a Steiner 
spread whose members are all S(2,n,g)s. □ 
The above theorem is useful to obtain regular LDPC 
codes through Theorems [T8l and [T9l and similar subdesgin 
deletion techniques based on Theorem [8] One can also 
modify Theorem [21] for the case when a GDD has dif- 
ferent group sizes by a similar argument. The existence 
of GDDs and their constructions have been extensively 
investigated in combinatorial design theory. For a com- 
prehensive list of known existence results on GDDs, we 
refer the reader to 14711 . 
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III. FINITE GEOMETRY CODES 

In this section, we demonstrate applications of our gen- 
eral designing methods by using combinatorial designs 
arising from finite geometries. 

The classical LDPC codes obtained from finite geome- 
tries are known to have remarkable error correction abil- 
ities. By using these codes, we generate infinitely many 
new high performance entanglement-assisted quantum 
LDPC codes having numerous Steiner spreads of vari- 
ous sizes. The various Steiner spreads in each code allow 
the code designer to flexibly fine-tune the parameters and 
error correction performance. 

This section is divided into three subsections. Subsec- 
tion [HO] studies entanglement-assisted quantum LDPC 
codes of girth six obtained from projective geometries. 
Codes based on affine geometries are investigated in Sub- 
section HTlB] In Subsection IIII CI we investigate slightly 
modified affine geometry codes, called Euclidean geom- 
etry codes. Classical LDPC codes based on these three 
kinds of finite geometries are called finite geometry LDPC 
codes or simply FG-LDPC codes. 

Many of the results presented in this section can also 
be seen as new results on classical finite geometry LDPC 
codes. In particular, properties of finite geometries have 
been independently studied in the combinatorial litera- 
ture, and hence many of the "known" results are new 
results in the field of LDPC codes. For the convenience 
of the reader, we summarize our results on fundamen- 
tal parameters of LDPC codes from finite geometries in 
Tables IXIVI and IXVI in Appendix [B] Lengths, dimen- 
sions, and minimum distances of the FG-LDPC codes 
with girth six from projective geometry PG{m, q), affine 
geometry AG(m,q), and Euclidean geometry EG{2,2 t ) 
are all determined. Specifically for EAQECCs based on 
FG-LDPC codes, we also determine the required amounts 
of entanglement for most cases. For a few cases, we give 
upper bounds on the required amount of entanglement. 



A. Projective geometry codes 

We begin with EAQECCs obtained from finite pro- 
jective geometries. The use of projective geometries for 
constructing EAQECCs first appeared in the work of 
Hsieh, Yen, and Hsu [l7j • This subsection illustrates how 
our combinatorial framework generalizes their method 
and determines fundamental parameters of quantum and 
classical LDPC codes obtained from PG(m, q). 

Points of the m-dimensional projective geometry 
PG(m, q) over ¥ q are the 1-dimensional subspaces 
of F™ +1 . The i-dimensional projective subspaces of 
PG(m,q) are the (i + l)-dimensional vector subspaces 
of F™ +1 . The points and lines of PG{m,q) form 

an S(2,q+ 1, q ^j -1 ), denoted by PGi(m,q), having 

L3. -0 (ff l) Win site o n A t-ot-»1 t r»o + i nn nnmKar 3. i — 



One can obtain two types of EAQECCs from projec- 
tive geometry designs: Type II (using a point-by-block 
incidence matrix) and Type I (using a block-by-point in- 
cidence matrix of the design) . Applying Proposition [3] 
to an incidence matrix of PGi(m,q), we obtain a Type 
II EAQECC. This type of EAQECC belongs to the high 
rate entanglement-assisted quantum LDPC codes given 
in Theorems [7] and [5] If we apply Proposition @] to 
a block-by-point incidence matrix, we obtain a Type I 
EAQECC. This kind of EAQECC belongs to the high 
redundancy entanglement-assisted quantum LDPC codes 
given in Theorem [§1 

The rank of an incidence matrix determines the dimen- 
sion of the corresponding FG-LDPC code, hence it is one 
of the key values in the quantum setting as well. Exact 
values for many sporadic examples have been computed 
in the fields of quantum and classical LDPC codes. The 
following two theorems give the exact rank for all projec- 
tive geometry designs. 



Theorem 22 (Hamada [38]) The rank o/PGi(m,2*) 
is given by 



rankPGi(m, 2*) = y>(m,2*) = 

i-li(s 3 +i,Sj) / , i\ / , o o- 

m + 1 \ (m+ 2sj + i — Sj — ii 
i I \ m 



e n e *' 

(s o ,si,...,s«)i=0 1=0 



where I = — 1, the sum is taken over all ordered sets 
(so, Si,..., St) with so = St, Sj 6 Z such that < Sj < 
m—1 andO < 2sj+i—Sj < ro+l for each j = 0, . . . ,t— 1, 
and 



2s j+1 - Sj 



blocks and replication number 
+ + 1. 
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L(s j+ i,Sj) 



We will use the notation (f(m, 2*) for the rank of 
PGi(m,q) when q is even, that is, q — 2*. When q 
is odd, the rank of PGi(m,q) is given by a formula of 
Frumkin and Yakir |5Q | . 

Theorem 23 (Frumkin and Yakir (50t ] ) Let q be odd 

and H an incidence matrix of the design PGi(m, q) with 

v = f-!- points. Then rank H = v — 1 = +^ . 

9-1 y <?-i 

Hence the exact dimensions of the corresponding FG- 
LDPC codes obtained from projective geometries can be 
calculated for all cases. 

The rank of PGi(m, 2*) was conjectured by Hamada 
[35l | to be the lowest rank among all Steiner 2-designs of 
the same order and block size. This has been confirmed 
in a number of cases, although in general the conjecture 
is still open. Thus we expect that the designs PGi(m, 2*) 
should provide codes with the best possible dimensions 
among all non-isomorphic S(2, 2* + 1, 2 { 2 t_\^ 1 )s. 

We will now examine the codes obtained from 
PGi(m,q) in detail. This subsection is divided into two 
portions based on the orientation of the incidence matrix. 
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1. Point-by-block (Type II) Projective geometry codes 

In this portion, we consider the EAQECCs correspond- 
ing to a point- by-block incidence matrix of PGi(m,q). 

We first consider the case q = 2* for some positive in- 
teger t. The following theorem gives an infinite family of 
entanglement-assisted quantum LDPC codes which con- 
sume only one initial ebit and have extremely large net 
rate. 

Theorem 24 For every pair of integers t > 1 and m > 
2 there exists an entanglement- assisted quantum LDPC 
codes with girth six whose parameters [[n, fc, d; c]] are 

(2*( m + 1 ) - l)(2 tm - 1) 

71 = 

(2 2t - 1)(2* - 1) ' 

(nt(m+l) _ lVO*m _ 1 "| 

d = 2 l + 2, and 
c= 1. 

To prove Theorem 1241 we first prove a new result on 
the distance of EAQECCs obtained from an incidence 
matrix of PGx(m, 2'). We use a special set of lines. A 
dual hyperoval TL is a set of q + 2 lines of PG i(2, q), such 
that each point of PG\{2 1 q) lies on either zero or two 
lines of H. Dual hyperovals exist if and only if q is even. 
An example is the set of projective lines with equations 

{A + pXi + !3 2 X 2 = : p e ¥ q } U {X x = 0} U {X 2 = 0}. 

Theorem 25 Let H be an incidence matrix of 
PGi(m, 2 ). TTie minimum distance of the classi- 
cal binary linear code with parity-check matrix H is 
2 f + 2. 

Proof. First, we note that coordinates of the code- 
words correspond to lines of the geometry, and a code- 
word corresponds to a set S of lines in PGi(m, 2*) such 
that every point is contained in an even number of lines 
of S. Assume that c is a non-zero codeword, and let 
supp(c) denote the support of c, that is, the set of in- 
dices of the nonzero coordinates of c. Since c ^ 0, the 
support of c contains at least one line £. Through each 
point of PG(m, 2*), there pass an even number of lines 
from supp(c). In particular, each of the 2* + 1 points on 
£ lies on at least one other line of supp(c), and all these 
lines are different as they have different intersections with 
£. Hence, there are at least 1 + (2* + 1) lines in supp(c), 
that is, minimum distance d is at least 2* + 2. Let 7r be a 
plane in PG(m, 2*) and S the set of the 2* + 2 lines of a 
dual hyperoval in n. Then S corresponds to a codeword 
of weight 2* + 2, hence d = 2* + 2. □ 



Proof of Theorem 1241 Let H be an incidence matrix of 
PGi(m, 2*). The rank of H is f{m, 2*) given by Theorem 
|2"2"1 The index of PG\(m ) 2 t ) is one. The replication 
number is odd. By Equation @ and TheoremO we have 
rank HH T = I. By Theor cm 1251 the minimum distance 
of the binary linear code with parity-check matrix H is 
2 f + 2. □ 
Next, we examine EAQECCs obtained from an in- 
cidence matrix of PGi(m,q) with q odd. This case 
also gives very high rate entanglement-assisted quantum 
LDPC codes. 

Lemma 26 Let H be an incidence matrix of PGx(2,q), 
q odd. Then the classical binary linear code defined by 
parity-check matrix H consists of only the zero vector 
and the all- one vector. 

Proof. This follows directly from Theorem |2"51 

A hyperbolic quadric Q is a substructure (P,C) of 
PGi(3,q) with (q + l) 2 points and 2(q + 1) lines, such 
that each point of V lies on exactly two lines of £ and 
every plane of PG(3, q) contains zero or two lines of C. 
Hyperbolic quadrics exist for every odd prime power q. 

Theorem 27 Let H be an incidence matrix of 
PG\{m,q), m > 3, q odd. Then the minimum 
distance of the classical binary linear code with a 
parity-check matrix H is 2(q + 1). 

Proof. Let n be a 3-dimensional subspace of PG(m, q) 
and (V, C) a hyperbolic quadric in n. The set of lines C 
determines a codeword of weight 2q-\-2, since each point 
of PG(m, q) is contained in zero or two lines of C. Hence 
minimum distance d is at least 2q + 2. 

We show that there are no codewords of weight smaller 
than 2q + 2. Assume that there exists a codeword c of 
weight smaller than 2q + 2, that is, supp(c) is a set of less 
than 2q + 2 lines of PG(m, q), such that each point lies on 
an even number of lines of supp(c). We will show that for 
any 2-dimensional subspace tt one has either |supp(c) n 
tt\ < 1 or |supp(c) (~1 7r > q + 2. 

First, let S = supp(c) D 7r = {£i, . , . ,£i}. For each 
j 6 {1, . . . , i}, each of the points on lj has to lie on at 
least one other line of supp(c), and at most i — 1 of them 
can lie on a line of S. Hence, at least q + 1 — (i — 1) 
of them are lines in supp(c) \ S and since they all have 
different intersections with tt, this yields i(q — i + 2) lines 
in supp(c) \ S. Together with the i lines of S, we have 

i(q - i + 2) + i < 2q + 2 

and solving this quadratic inequality for i gives us that 
either i>q + lori<2. Since i is an integer, hence 
i > q + 2 or i < 1. 

Now, let I be any line of supp(c). Each point of £ must 
lie on at least one other line, hence there certainly exist 
planes n with i > 2, and we have i > q + 2. Let n be such 
a plane. We will now show that all lines of supp(c) are 
contained in n. Assume the contrary, that there exists 
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a line t' € supp(c) \ S. Through each of the points on 
£' \ 7r, we need at least one other line of supp(c) which is 
not contained in tt. Since there are at least q points on 
£' \ 7T, one has 

|supp(c)| = |5| + |supp(c)\S , | > (q + 2) + (l + q) > 2q + 2, 

a contradiction. Hence, £' does not exist and supp(c) 
is contained within a single plane tt. However, tt is a 
PGi(2, q) and by Lemma [26] we need q 2 + q+l>2q + 2 
lines in this contradiction. Hence, there are no 

codewords of weight less than 2q + 2. □ 
We now give another infinite family of Type II 
entanglement-assisted quantum LDPC codes. 

Theorem 28 Let q be an odd prime power. Then for ev- 
ery integer m > 3 there exists an entanglement-assisted 
quantum LDPC code with girth six whose parameters 
[[n, k, d\ c]] are 

fq m+l _ ^fqm _ X ) 

ft 

(<Z 2 -l)(g-l) ' 



_ [q m + l - l)(q m - 1) _ n q m+1 -q 

( q *-l)(q-l) q-1 +C ' 

d = 2q + 2, and 

{1 when m is odd, 
q g _ 1 ~ <? when to is even. 

Proof. This follows directly from Proposition [3] and 
Theorems HESl and EH □ 
Therefore in the case where m is odd, we have another 
infinite class of EAQECCs which consume only one ebit. 
If 777 is even, we obtain infinitely many high rate codes 
which consume reasonable numbers of ebit. Tables U and 
llTl give a sample of the parameters of the Type II codes 
obtained from PGi(m,q) with q even and q odd respec- 
tively. 



TABLE I 




Sample parameters 


of Type 


II [[n, k, o 


!;c]] 


EAQECCs 


obtained from PGi(m 


,q), q even. 






m 


q 


ii 


k 


d 


c 


3 


2 


35 


14 


4 


1 


4 


2 


155 


104 


4 


1 


5 


2 


651 


538 


4 


1 


6 


2 


2667 


2428 


4 


1 


3 


4 


357 


236 


6 


1 


4 


4 


5795 


5204 


6 


1 


2 


8 


73 


18 


10 


1 


3 


8 


4745 


3944 


10 


1 



In the reminder of this portion, we examine Steiner 
spreads of projective geometry designs. These substruc- 
tures can be used in Theorems [T71 [T51 and Ql|] and their 



TABLE II. Sample parameters of Type II [[n, k, d; c]] 
EAQECCs obtained from PGi(m,q), q odd. 



771 


Q 


n 


k 


d 


c 


3 


3 


130 


53 


8 


1 


3 


5 


806 


497 


12 


1 


3 


7 


2850 


2053 


16 


1 


4 


3 


1210 


1090 


8 


120 



analogous techniques based on Theorem [5] to fine-turn 
the rates and distances of the EAQECCs. 

An s-spread of PG(m 7 q) is a set of s-dimensional pro- 
jective subspaces which partition the points of the ge- 
ometry. In other words, an s-spread consists of a set of 
(,s + l)-dimensional vector subspaces of F™ +1 which con- 
tain every nonzero vector exactly once. It is known that 
PG(m, q) admits an s-spread if and only if s + 1 divides 
777 + 1 (see [5lf] and H p. 29]). 

Take PGi(m,q) and suppose s > 2 is chosen so that 
s + 1 divides m + 1. Then an s-spread of PG(to, q) exists. 
Each s-dimcnsional subspace in the spread contains an 
isomorphic copy of PGi(s, q), and hence this forms a 
Steiner spread. Note that the blocks of PGi(s,q) have 
size q + 1 and are also blocks of PGi(m, q). Therefore we 
have the following result. 

Theorem 29 Let s, m > 1 be positive integers such that 
s + 1 divides m + 1. Then PGi(m,q) contains q qS +i Z\ 
disjoint copies of PG\(s, q) whose point sets partition the 
point of PGi(m,q). 

Thus, we can find a set of disjoint subdesigns which par- 
tition the points of PG\ (m, q) whenever m + 1 has a non- 
trivial factor. Naturally, we may further sub-divide each 
subdesign of dimension s into smaller subdesigns, based 
on the nontrivial factors of s + 1. Hence, the S(2, /j, v)s 
from PGi (m, q) are very flexible in that they have Steiner 
spreads of various sizes. 

In general, the length, dimension, required ebits, and 
rate each change gradually as we delete subdesigns in 
a Steiner spread. The minimum distance and rank are 
either remain the same or improve slightly. Table IIIII 
lists the example parameters of EAQECCs created by 
deleting subdesigns from PGi(5,2). The first and last 
rows correspond to regular LDPC codes. 



2. Block-by-point (Type I) Projective geometry codes 

Next we consider EAQECCs obtained via Theo- 
rem [9] by using the block-by-point incidence matrix of 
PG\(m,q). The codes obtained in this manner corre- 
spond to the classical Type I LDPC codes. As in the 
classical setting, Type I entanglement-assisted quantum 
regular LDPC codes can correct many quantum errors. 
Because an incidence matrix of PG\(m,q) for q odd is 
almost full rank, the corresponding Type I code is not of 
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TABLE III. Summary of Type II codes obtained by deleting 
a Steiner spread of subdesigns isomorphic to PGi(2,2) from 
PGi(5,2). 



Subs a 


71 


ran W T~f 

L Ctllxv J. J. 


k 






Rate 


n 
u 


Uj ± 


^7 

O i 


ooo 


A 
*± 


1 
1 




1 


644 


57 


532 


4 


2 


0.8370 


2 


637 


57 


526 


4 


3 


0.8477 


3 


630 


57 


520 


4 


4 


0.8587 


4 


623 


57 


514 


4 


5 


0.8700 


5 


616 


57 


508 


4 


6 


0.8815 


6 


609 


57 


502 


4 


7 


0.8933 


7 


602 


57 


496 


4 


8 


0.9053 


8 


595 


57 


490 


4 


9 


0.9176 


9 


588 


57 


482 


4 


8 


0.9269 



a This column denotes the number of subdesigns removed. 

much interest. Hence, in this portion we always assume 
that q = 2* for some positive integer t. 

Theorem 30 For every pair of integers t > 1 and m > 
2 there exists an entanglement- assisted quantum LDPC 
code with girth six whose parameters [[n, k, d] c]] are 

2 t(m+l) _ 1 



2 t(m+l) _ x 

fc = y _ 1 2<p(™,2 t ) + c, 

d= (2* + 2)2* (m - 2) , and 
c < ^{m,2 l ). 

Proof. Let H T be a block-by-point incidence matrix 
of PGi(m,2*). Then we have rank# T # < rankiJ = 
<p{m, 2*), where <p(m, 2*) is given by Theorem |2"2"1 By a 
result of Calkin, Key, and de Resmini [45| , the minimum 
distance of the binary linear code with parity-check ma- 
trix H T is (2* + 2)2* (m ~ 2) . Applying Propositiong]proves 
the assertion. □ 

Note that here the distance grows exponentially as the 
dimension of the geometry increases. When m = 2, the 
EAQECCs are based on projective planes. As shown in 
Subsection III B[ the EAQECC obtained from a Desar- 
guesian projective plane of order 2* consumes only one 
initial ebit. Basing on Hamada's conjecture, we expect 
that in general the EAQECCs given in Theorem [501 con- 
sume relatively small numbers of ebits. 

It is not clear from the formula for </j(to, 2') whether 
the net rate of a Type I EAQECC based on PGi(m, 2*) 
is positive. In order to produce useful catalytic quantum 
codes, it is important to understand when the net rate is 
positive. 

Proposition 31 Let H be an incidence matrix of 
PGi(2,2 t ). Then for all t > 2 the EAQECC obtained 
from H T has a positive net rate. 



Proof. By Hamada's formula, we have ranki7 = 3*+l. 
The number of points in PGi(2, 2*) is 2 2t + 2* + 1. □ 
For m > 3, we note that as q increases, ranki? grows 
at a slower rate than the code length. Thus we may 
expect that, for q large when compared to m, the net rate 
will eventually become positive. For example, one can 
check that the net rate of the Type I EAQECC obtained 
from PGi(3,2*) is positive for t > 7. Table [TV] gives 
sample parameters of the Type I codes obtained from 
PGi(m,2*). 



TABLE IV. Sample parameters of Type I [[n, k, d; c]] 
EAQECCs obtained from PGi(m,q), q even. 



in 


q 


n 


k 


d 


c 


2 


4 


21 


2 


6 


1 


2 


8 


73 


18 


10 


1 


2 


16 


273 


110 


18 


1 


2 


32 


1057 


570 


34 


1 



B. Affine geometry codes 

In this subsection, we will study the EAQECCs ob- 
tained from affine geometry designs. 

The affine geometry AG(m, q) of dimension m over ¥ q 
is a finite geometry whose points are the vectors in F™. 
The i-dimcnsional affine subspaces (or i- flats) are the 
i-dimensional vector subspaces of F™ 1 and their cosets. 
Thus AG(m, q) has a natural parallelism. 

The points and lines (that is, 1-flats) of an affine 
geometry form an S(2,q,q m ), denoted by AG\{m, q). 
The design has q m ~ x q q I^ blocks and replication num- 
ber S^fii = q ™-i + q m - 2 + h q + 1. 

We note that in many papers concerning LDPC 
codes, the term "Euclidean geometry" and the notation 
EG(m, q) are used for affine geometries. Most of the 
codes studied in relation to Euclidean geometries does 
not use the zero vector, and hence they do not gener- 
ally correspond to 5(2, /i, v)s. Because the term "afhnc 
geometry" is standard in the recent research on finite ge- 
ometry in mathematics, we use the notation AGi(m,q) 
when we take all points and lines to form an incidence 
matrix. The incidence structure obtained by discarding 
the zero vector and the lines containing the zero vector 
from AGi(m, q) will be denoted by EGi{m, q), which we 
will study in Subsection IIII CI Because many of the clas- 
sical FG-LDPC codes obtained from affine geometries are 
based on EG\{m, q), they are generally not the same as 
the affine geometry codes presented in this section. 

As with projective geometry designs, Propositions [3J 
and H] give Type II and Type I affine geometry codes re- 
spectively. It is notable that the classical ingredients of 
these codes are quasi-cyclic LDPC codes similar to other 
FG-LDPC codes because the elementary abeli an grou p 
acts transitively on the points of AG\(m, q) (see |18l.l32j|). 
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The rank of an affine geometry design AGi(m, 2') is di- 
rectly related to tp given in Theorem |2"2"1 

Theorem 32 (Hamada |35l |) The rank of the affine 
geometry design AGi(to,2') is given by 

rank AG x (m, 2*) = ¥?( m > 2 *) - p(m - 1,2*). 

If q is odd, the rank of AG\(m, q) over F2 is full. 

Theorem 33 (Yakir [53]) Let H be an incidence ma- 
trix of the design AG\(m,q) with q odd. Then ranki? = 
q m . 

Thus the dimensions of the corresponding FG-LDPC 
codes can be easily calculated. 

As in the case of projective designs, Hamada conjec- 
tured that the rank of AG\(m, 2*) is minimum among all 
Steiner 2-designs of the same order and block size. Thus, 
affine geometry designs with q even may be expected to 
give codes with the best possible dimensions among all 
non-isomorphic S(2, 2*, 2 tm )s. 

We divide this subsection into two portions. In 
the first portion we examine high rate Type II 
entanglement-assisted quantum LDPC codes obtained 
from AGi(m,q). Then in the next portion we present 
Type I entanglement-assisted quantum LDPC codes 
based on AGi(m,q), which effectively exploit the redun- 
dancy to give excellent error correction performance. 

1. Point-by-block (Type II) Affine geometry codes 

The geometric structure of affine geometry has often 
been studied independently in various fields. The spe- 
cial substructure we need to give distances has been in- 
vestigated in connection with the disk failure resilience 
ability of a class of redundant arrays of independent disks 
(RAID). Here we present a known result on RAID related 
to our codes in coding theoretic terminology. 

Theorem 34 (Miiller and Jimbo [54]) Let H be an 

incidence matrix of AG\{rn,q). The minimum distance 
of the classical binary linear code having H as a parity- 
check matrix is q + 1 if q is even, and 2q otherwise. 

The following two theorems give infinite families of 
EAQECCs which consume only one initial ebit and have 
very large net rate. 

Theorem 35 For every pair of integers t > 1 and m > 
2 there exists an entanglement- assisted quantum LDPC 
code with girth six whose parameters [[n, k, d; c]] are 



cytrn i 

k = 2 t (" l - 1 )^ r - r - 2( ¥ >(m, 2 4 ) - <p(m - 1,2')) + 1, 

d = 2* + 1, and 
c= 1. 



Proof. Let H be an incidence matrix of AGi(m, 2'). 
By Theorem 132] we have rank_ff = (f[m, 2 ) — tp(m — 
1,2*). The index of the design AGi(m,2 t ) is one. Its 
replication number is always odd. Thus, by Theorem [7] 
we have rank HH T = 1 . Applying Proposition [3] and 
Theorem l34l completes the proof. □ 

Theorem 36 Let q be an odd prime power. Then for ev- 
ery integer to > 2 there exists an entanglement- assisted 
quantum LDPC code with girth six whose parameters 
[[n, k, d; c]] are 



in, _ 

k = q m -^ 2q m + c, 

9-1 

d = 2q, and 

{1 when m is odd, 
q m — 1 when m is even. 

Proof. Let H be an incidence matrix of AG\ (to, q) with 
q odd. By Theorem 1331 we have ranki? = q m . The in- 
dex of the design AG\ (to, q) is one. Its replication num- 
ber r is a sum of to terms, each being an odd number. 
Thus r is odd only when to is odd. By Theorem [71 we 
have iankHH T = 1 for to odd. If to is even, we have 
rank HH T = q rn — 1 from Theorem [H Applying Propo- 
sition |3] and Theorem [M] proves the assertion. □ 
Theorem [36] gives an infinite family of high rate 
entanglement-assisted quantum LDPC codes which ex- 
ploit reasonable amounts of entanglement as well. Tables 
IVl and I VII give a sample of the parameters of the Type 
II codes obtained from AG\(m, q) with q even and q odd 
respectively. 



TABLE V. Sample parameters of Type II [[n, k, d; c]] 
EAQECCs obtained from AGi(m,q), q even. 



rn 


q 


n 


k 


d 


c 


3 


2 


28 


15 


3 


1 


4 


2 


120 


91 


3 


1 


5 


2 


496 


435 


3 


1 


6 


2 


2016 


1891 


3 


1 


2 


4 


20 


3 


5 


1 


3 


4 


336 


235 


5 


1 


4 


4 


5440 


4971 


5 


1 


2 


8 


72 


19 


9 


1 


3 


8 


4672 


3927 


9 


1 



Next we show that affine geometry designs have nu- 
merous subdesigns and Steiner spreads, which make it 
possible to fine-tune the parameters and error correction 
performance of the corresponding EAQECCs. 
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TABLE VI. Sample parameters of Type II [[n, k,d;c]] 
EAQECCs obtained from AGi (m, q) , q odd. 



m 


q 


n 


k 


d 


c 


3 


3 


117 


64 


6 


1 


3 


5 


775 


526 


10 


1 


3 


7 


2793 


2108 


14 


1 


5 


3 


9801 


9316 


6 


1 


4 


3 


1080 


998 


6 


80 



Theorem 37 If m > 3, the points of AGi(m,q) can be 
partitioned into q disjoint subsets of size q m ~ x , being the 
point sets of subdesigns isomorphic to AG\(m — l,q). 

Proof. Take a parallel class {Hi,...,H q } of q hy- 
perplanes of AG(m,q). Let the point set of Hj be 
Vj. Clearly \Jj =1 Vj — V, and the set of all blocks of 
AG\ (m, q) which are contained entirely in Hj form a sub- 
design isomorphic to AG\{m — 1, q). □ 
Theorem 1371 can be applied recursively to create addi- 
tional disjoint subdesigns of smaller dimension, giving a 
variety of EAQECCs via Theorems El US andQH Sim- 
ilar subdesign deletion techniques based on Theorem [5] 
further give infinitely many new high rate EAQECCs. 
Table IViTl lists the parameters of the EAQECCs created 
by spread deletion from AC?i(3,4). 



TABLE VII. Summary of Type II codes obtained by deleting 
a Steiner spread of subdesigns isomorphic to ^4Gi(2,4) from 
AGi(3,4). 



Subs^ 


n 


rank H 


k 


d 


c 


Rate 





336 


51 


235 


5 


1 


0.6994 


1 


316 


51 


216 


5 


2 


0.7468 


2 


296 


51 


197 


5 


3 


0.8007 


3 


276 


51 


178 


5 


4 


0.8623 


4 


256 


51 


158 


6 


4 


0.9297 



a This column denotes the number of subdesigns removed. 



2. Block-by-point (Type I) Affine geometry codes 

Next we consider the EAQECC obtained from a block- 
by-point incidence matrix of AG\(m, q). Because in- 
cidence matrices of AG\(m, q) with q odd are of full 
rank, here we always assume q = 2* to obtain interest- 
ing codes. The entanglement-assisted quantum LDPC 
codes presented in this section effectively exploit redun- 
dancy. The excellent error correction performance will 
be demonstrated in simulations in Section IIV1 

Theorem 38 (Calkin, Key, and de Resmini [45]) 

Let H be a block-by-point incidence matrix of AGi(m, 2*). 
Then the minimum distance of the classical binary lin- 
ear code for which H is a parity-check matrix is 
(2* + 2)2 t( - m - 2 K 



Theorem 39 For every pair of integers t > 1 and to > 
3 there exists an entanglement- assisted quantum LDPC 
code with girth six whose parameters [[n, fc, d; c]] are 



k = 2 tm - 2((p(m, 2*) - ip(m - 1, 2*)) + c, 
d = (2* + 2)2*(" l ~ 2 \ and 

c < vim^*) - ip{m - 1,2*). 

Proof. Let H T be a block-by-point incidence matrix 
of AGi(m,2*). By Theorem [32 we have mnk H T H < 
rankff = y>(m,2*) - <p(m - 1,2*). By Theorem |3"51 the 
minimum distance of the binary linear code with a parity- 
check matrix H is (2* + 2)2*( m ~ 2 ). The assertion follows 
from Proposition |4j □ 
It is worth mentioning that here the distance grows ex- 
ponentially with linear increase of the geometry dimen- 
sion to. Because the rank of AG\(m, 2*) is conjectured 
to be the smallest possible among all non-isomorphic 
5(2,2*,2 tm )s, we expect that the EAQECCs obtained 
from these affine geometry designs consume the smallest 
possible numbers of ebits attainable by this method with 
S'(2,2 t ,2 tm )s. 

When m = 2, we can easily determine the required 
amount of entanglement. 

Theorem 40 For every positive integer t there exists an 
entanglement- assisted quantum LDPC code with girth six 
whose parameters [[n, k, d; c\] are 

n = 4', 
fc = 4* + 2*-2-3*, 
d = 2* + 2, and 
c = 2*. 

Proof. Let H T be a block-by-point incidence matrix of 
AG X (2,2*). We first prove that rankiJ T iJ = 2*. Two 
lines of an affine plane are either parallel or intersect in 
exactly one point. There are 2* + l parallel classes of lines, 
each containing exactly 2* lines, and each line contains 2* 
points. Because 2* is even, it is always possible to reorder 
the rows of H T such that H T H is a block matrix of the 
following form: 

" J J' 



J J ■■■ 

where J is the 2* x 2* all-one matrix. Hence, we have 
rank H T H = 2*. By Theorem W\ we have rank H = 3 . 
Applying Proposition [4] and Theorem [38] completes the 
proof. □ 
Table IVIIII gives sample parameters of the Type I 
EAQECCs obtained from AG\[m, 2'). 
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TABLE VIII. Sample parameters of Type I [[n, k, d; c]] 
EAQECCs obtained from AGi (m, q) , q even. 



Ill 


q 


n 


k 


d 


c 


2 


8 


64 


18 


10 


8 


2 


16 


256 


110 


18 


16 


2 


32 


1024 


570 


34 


32 


2 


64 


4096 


2702 


66 


64 



C. Euclidean geometry codes 

In this final subsection concerning finite geometry 
EAQECCs, we will examine Euclidean geometry codes. 

Given a prime power q and integer m > 2, we define an 
incidence structure EG\ (to, q) having as points all points 
of AGi (to, q) except the zero vector, and having as blocks 
(or lines) all lines of AG(m, q) except those lines contain- 
ing the zero vector. The lines which are excluded from 
AGi(m, q) to form EGi(m, q) consist of all multiples of a 
single nonzero vector. Thus, EG\(m, q) has q m — 1 points 
and (q m ~ x — l) q q Z\ lines. Each line contains q points, 

and each point appears in q q ~i ~~1 = q m ~ 1 +q m ~ 2 +- ■ -+q 
lines. Thus, EGi(m, q) yields regular LDPC codes. Each 
pair of points appears in at most one line. Hence, 
EG\{m, q) is a partial Steiner 2-design. Its Tanner graph 
does not contain 4-cycles. 

Applying Proposition U to a line-by-point incidence 
matrix of EG\(m,q) gives a Type I EAQECC. If q is 
even, the distance is bounded from below by the BCH 
bound. 

Theorem 41 (Kou, Lin, and Fossorier |l8j ]) Let H 

be a line-by-point incidence matrix of EG±(m, 2*). Then 
the minimum distance d of the classical binary linear code 
having H as a parity-check matrix satisfies d > \ t . 
Equality holds if m = 2 . 

We use the following theorem to give the dimensions 
of FG-LDPC codes obtained from EGi(m, 2') and their 
entanglement-assisted quantum counterparts. 

Theorem 42 (Hamada |35l ] ) The rank of the inci- 
dence structure -EGi(m, 2'), t > 1, is given by 

rank EGi(m, 2*) = ¥>(m,2*) - tp(m - 1,2*) - 1. 

Theorem 43 For every pair of integers t > 1 and m > 
2 there exists an entanglement- assisted quantum LDPC 
code with girth six whose parameters [[n, k, d; c]] are 

n = 2 tm - 1, 
k = 2 tm - 2(<p(m, 2*) - <fi(m - 1, 2*)) + 1 + c, 

d > — : , and 

~ 2* - 1 ' 

c < ip(m, 2*) - ip(m - 1, 2') - 1. 



Proof. Let H T be a line-by-point incidence matrix of 
EG 1 (m,2 t ). By Theorem SI we have rank H T H < 
ranki7 = (f(m, 2*) — <p{m — 1, 2') — 1. Applying Propo- 
sition S] and Theorem WW completes the proof. □ 
A simple observation gives exact values of all the pa- 
rameters of the Type I codes based on EG\{2, 2'). 

Theorem 44 For every positive integer t there exists an 
entanglement- assisted quantum LDPC code with girth six 
whose parameters [[n, k, d; c)] are 

n = 4* - 1, 
fc = 4* + 2 t -2-3 t + l, 
d = 2* + 1, and 
c = 2 f . 

Proof. Let H T be a line-by-point incidence matrix of 
£Gi(2,2*). An incidence matrix of EGi{2,2 t ) is ob- 
tained by removing one row and one column from each 
block from that of AGi(2,2 t ). By following the argu- 
ment in Theorem SOI h is straightforward to see that 
rank H T H = 2*. By Theorem we have ranki? = 
(f(m, 2') - (p(m - 1, 2') - 1 = 3* - 1. Theorem 04] and 
Proposition SI prove the assertion. □ 
Table ITXl gives a sample of the parameters of the Type 
I codes obtained from EG\{2, 2*). 



TABLE IX. Sample parameters of Type I [[n, k, d; c]] 
EAQECCs obtained from EG\(2,q), q even. 



m 


q 


n 


k 


d 


c 


2 


8 


63 


19 


9 


8 


2 


16 


255 


111 


17 


16 


2 


32 


1023 


571 


33 


32 



As with S(2, fi, v)s, the incidence structure EGi(m, q) 
can also generate a high rate LDPC code with girth six. 
Applying Proposition [3] to incidence matrices, we obtain 
Type II EAQECCs. Here we investigate their parame- 
ters. 

Theorem 45 The minimum distance of a Type II 
EAQECC based on EG\{m,q) is q+ 1 if q is even, and 
2q if q is odd and m > 2 . 

Proof. Consider any set of linearly dependent columns 
in an incidence matrix of EGi(m,q). The same 
columns appear in the corresponding incidence matrix 
of AG\{m, q), but with a single zero coordinate added. 
These columns are still dependent in AG\{m,q). Hence 
the minimum distance is upper bounded by Theorem [34] 
Thus we need only to show lower bounds. 

We begin with q even. If q = m = 2, we can check 
by hand that the minimum distance is three. Henceforth 
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assume that q > 2 or m > 2. Because the minimum 
distance of the code obtained from AGi(m,q) is q + 1, 
there exists a set S of q + 1 linearly dependent columns 
of an incidence matrix of AGi(m,q), corresponding to a 
set V of q + 1 blocks of AG\ (to, q). Let P be the multiset 
of points appearing in the blocks of T>. As each block 
of T> has q points, \P\ = q(q + 1). However, because 
the columns of S are dependent over F2, each point in P 
must appear with multiplicity two or more. Hence, the 
number of distinct points in P is at most q ( q + 1 ' 1 < q m — l 
except for q = m = 2. Therefore there is a nonzero 
point p of ^4G(m, q) which does not appear in P. Let 
V = {B - p : B eV], that is, we shift each block of V 
by p. Each new block corresponds to a coset of a linear 
space. Because p P, no element of 2?' contains the zero 
vector, and so the elements of V are lines of EGi(m, q). 
Thus V is a linearly dependent set in EG\{m,q) of size 
(7 + 1. Therefore in all cases, the minimum distance of 
Type II EAQECC based on EG\{m,q), q even, is q + 1. 
A similar argument proves the case when q is odd and 
rn ^ 2. □ 

Theorem 46 For every pair of integers t > 1 and to > 
2 i/iere exists an entanglement- assisted quantum LDPC 
code im'f/i girth six whose parameters [[n, fc, 0!; c]] are 



(2 



t(m-l) 



1) 



1 



2* - 1 



(2 



t(m-l) 



2* - 1 



2rank J BGi(TO,2 t ) + c, 



2* + 1, ana 



2* - 1 

where ranki?Gi(TO, 2*) = (p(m, 2*" 1 



1. 



') - p(m - 1,2*) 

Proof. Let H be an incidence matrix of EG\(m : 2*). 
Because 7J is obtained from an incidence matrix of 
AG 1 (to, 2') by deleting the row representing the zero vec- 
tor and the columns that represent the lines containing 
the zero vector, it is easy to see that the rows and columns 
of HH T can be reordered such that the matrix is of the 
form: 



HH 1 



J 
J 



J J 



J 
J 







where J is the (2* - 1) x (2* 
2 tm _ i is oddj rank hh t = 

sition [3] and Theorems! 



- 1) all-one matrix. Because 
1 1 1. Applying Propo- 



2*-l 

] and |3H completes the proof. □ 
Tables IXl gives sample parameters for the Type II codes 
obtained from EG\(m, 2*). 

For the case q odd, Hamada [35j conjectured that an 
incidence matrix of EG\{rn,q) is of full rank. As shown 
in Table IXll the conjecture is true for small to and q. 



TABLE X. Sample parameters of Type 
EAQECCs obtained from EGi(m,q), q even. 



II \\n, k, d: 



m 


q 


n 


k 


d 


c 


3 


2 


21 


15 


3 


6 


4 


2 


105 


91 


3 


14 


5 


2 


465 


434 


3 


30 


6 


2 


1953 


1891 


3 


62 


3 


4 


315 


235 


5 


20 


4 


4 


5355 


4971 


5 


84 


2 


8 


63 


19 


9 


8 


3 


8 


4599 


3927 


9 


72 


TABLE 


XI. 


Sample parameters of Type 


II 


[n,k,d;c]} 


EAQECCs obtained from EGi( 


m,q), q odd. 






m 


q 


n 


k 


d 


c 


3 


3 


104 


64 


6 


12 


4 


3 


1040 


960 


6 


80 


5 


3 


9680 


9316 


6 


120 


3 


5 


744 


526 


10 


30 


3 


7 


2736 


2108 


14 


56 



IV. PERFORMANCE 

In this section, we present simulation results for 
EAQECC codes constructed in the previous sections. As 
in the related works [l(| H3] , we performed simulations 
over the depolarizing channel. In this model, each error 
(X, Y, and Z) occurs independently in each qubit with 
equal probability f m . For a given CSS type EAQECC, we 
performed each decoding in two separate decoding steps, 
each using the sum-product algorithm. The shared ebits, 
which do not pass through the noisy channel, are assumed 
to be error-free. Our simulation results are reported in 
terms of the block error rate (BLER) . 

We first examine codes obtained from a block-by-point 
incidence matrix. Figure [T] shows the performance of sev- 
eral such codes based on projective and affine geome- 
try designs. As shown in Section iHll these codes have 
very large distances for sparse-graph codes while avoid- 
ing short cycles. As expected, these codes perform excel- 
lently at relatively high f m . 

To illustrate how well these codes perform, we compare 
one of our Type I LDPC codes with previously known 
entanglement-assisted quantum LDPC codes with best 
BLERs. 

Theorem [40] gives a new EAQECC with parame- 
ters [[256,110, 18; 16]] obtained from AGi(2,16). The 
[[255, 111, 17; 16]] EAQECC in the work of Hsieh, Yen, 
and Hsu [13] used -EGi(2,16) outperformed all previ- 
ously known quantum codes of similar rate in simula- 
tions over the depolarizing channel. Their code based on 
PGi(2, 16), which also performed very well, has param- 
eters [[273,110, 18; 1]]. Exactly the same EAQECCs as 
these two can be constructed using Theorems |44] and \W\ 
in our framework without relying on computers to calcu- 
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FIG. 1. (Color online) Performance of Type I EAQECCs 



late their parameters. 

These three EAQECCs based on finite geometries have 
similar geometrical structures, and they behave quite 
similarly in simulations. Performance of the AGi(2, 16) 
and PGi(2,16) codes is directly compared in Figured] 
The BLER of the £Gi(2,16) code, which is slightly 
worse than that of our AGi(2,16) code, is plotted in 
Figure [2] to compare the three with EAQECCs having 
different parameters. As shown in the figures, our new 
[[256, 110, 18; 16]] EAQECC obtained from AGi(2,16) 
shows a better BLER than the other two. The BLERs of 
AGi(2,16), £Gi(2,16), and PGi(2,16) codes at f m = 
0.02 are 1.0 x 10~ 4 , 1.6 x 10~ 4 , and 3.8 x 10~ 4 respectively. 

Entanglement-assisted quantum quasi-cyclic LDPC 
codes proposed by Hsieh, Brun, and Devetak in [l6[ 
have also shown excellent BLERs. In simulations their 
[[128, 58, 6; 18]] EAQECCs, called EX1 and EX2, outper- 
formed the previously known best quantum LDPC codes 
at a similar rate about 0.316. The net rate of EX1 and 
EX2 is « 0.312. Our [[256, 110, 18; 16]] EAQECC 

obtained from AGi(2,16) has net rate "j?" 16 w 0.367, 
which is higher than that of EX1 and EX2. Their simu- 
lation results and our independent simulation results for 
EX1 and EX2 showed that their BLERs at f m = 0.02 
are higher than 1.1 x 10~ 2 while our AG\(2 : 16) code has 
BLER about 1.0 x 10~ 4 at the same f m , which is bet- 
ter than EX1 and EX2 by two orders of magnitude. Our 
EAQECC also requires a smaller amount of entanglement 



than EX1 and EX2. 

Our results here confirm the close linkage between 
EAQECCs and classical error-correcting codes: good 
performance in the classical setting translates directly 
into good performance from the corresponding quantum 
codes. 
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FIG. 2. (Color online) Performance of Type II EAQECCs 

We next examine codes obtained from a point-by-block 
incidence matrix. These codes are capable of achieving 
extremely high rates even at moderate block lengths. 

Figures [2] and [3] show the performance of several Type 
II codes based on finite geometries. The Type II code 
from PGi(3,3) is shown in both figures to serve as a 
point of reference between the two figures. Figure 2] gives 
the block error rates for several codes with high rates in- 
cluding [[301, 216, 6; 1]] and [[1080, 998, 6; 80]] codes from 
cyclic 5-sparse STSs of order 43 and 81 respectively. The 
incidence matrices of these two Steiner triple systems are 
constructed from the list of base blocks in [55j . Note 
that the cyclic automorphisms and sparse configurations 
immediately give the dimensions and distances of the 
EAQECCs obtained from the cyclic 5-sparse STSs (see 
[37l [56[ ) . Table IXIII lists the rates of selected finite ge- 
ometry codes shown in figures. 

As in the classical setting, our codes obtained from 
point-by-block incidence matrices have waterfall regions 
at low f m and transmit at extremely high rates. This di- 
rect correlation in performance between the classical and 
quantum settings can also be seen when codes require 



17 



OA 



0.01 



0.001 



0.0001 



0.00001 




0.001 0.003 0.005 0.007 

Crossover probability f m on the depolarizing channel 

FIG. 3. (Color online) Performance of Type II EAQECCs 



OA 



0.01 



0.001 



0.0001 









AG,(3,7) 
PG,(3,7) 
PG,(4,3) 




































5-sparse STS(43) 
















5-sparsc STSI81) 



























0.0005 0.0015 0.0025 0.0035 

Crossover probability f m on the depolarizing channel 



FIG. 4. (Color online) Performance of high-rate Type II 
EAQECCs 



TABLE XII. Rates of EAQECCs obtained from finite geome- 



tries. 


Type 


Geometry 


m 


q 


Rate 


II 


PG 


4 


3 


0.9008 


II 


PG 


3 


7 


0.7203 


II 


PG 


3 


5 


0.6166 


II 


PG 


3 


3 


0.4076 


II 


AG 


3 


7 


0.7547 


II 


AG 


3 


5 


0.6787 


II 


AG 


3 


3 


0.5470 


II 


AG 


2 


8 


0.2638 


II 


EG 


2 


16 


0.4352 


II 


EG 


2 


8 


0.3015 


I 


PG 


2 


32 


0.5392 


I 


PG 


2 


16 


0.4029 


I 


PG 


2 


8 


0.2465 


I 


AG 


2 


32 


0.5566 


I 


AG 


2 


16 


0.4296 


I 


AG 


2 


8 


0.2812 


only one 


ebit. It may 


be worth mentioning 


that chang- 



ing geometries or choosing a non-geometric 5(2, fi, v) can 
give slightly different BLER curves. It would be in- 
teresting to investigate theoretical methods for finding 
5(2, n, v)s with desirable performance curves in given sit- 
uations. 

Finally, we compare EAQECCs obtained by removing 
subdesigns from the parent design. Here we test a subde- 



sign deletion technique where each deletion step increases 
the required amount of entanglement to a slightly larger 
degree than the examples we gave in Section IIIII Each 
code in Figure[5]is constructed from a Type II code based 
on AG i (3,3). Fundamental parameters of these codes are 
shown in Table IXIIII The original code is also shown for 



TABLE XIII. Summary of Type II EAQECCs obtained by 
deleting subdesigns from AGi(3, 3). 



Subs a 


n 


rank H 


k 


d 


c 


Rate 





117 


27 


64 


6 


1 


0.5470 


1 


105 


27 


60 


6 


9 


0.5714 


2 


93 


26 


58 


6 


17 


0.6236 


3 


81 


25 


56 


6 


25 


0.6913 



a This column denotes the number of subdesigns removed. 



reference. The code labeled "one sub" has had a single 
subdesign isomorphic to AG i(2, 3) removed. The code 
labeled "3 subs" has had a Steiner spread removed. This 
last code is a regular LDPC code. As can be seen from 
their BLERs, removing subdesigns has improved the er- 
ror correction performance while increasing the rate and 
maintaining many of the essential properties. 

Because removing subdesigns can increase the required 
amount of entanglement in a flexible manner, one can 
generate an EAQECC which effectively exploits pre- 
existing entanglement. For example, a high net rate 
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FIG. 5. (Color online) Performance of EAQECCs obtained 
by deleting subdesigns from AGi(3, 3). 



code consuming only one ebit can turn into a heavily 
entanglement-assisted code to achieve better error cor- 
rection performance at the same f m . As illustrated in 
Table [XTTTl a [[117, 64, 6; 1]] code with a regular parity- 
check matrix becomes a [[81, 56, 6; 25]] code with a regu- 
lar parity-check matrix through gradual steps. 

One can also fine-tune parameters and improve er- 
ror correction performance while almost keeping the ex- 
tremely low required amount of entanglement by apply- 
ing Theorems [T7] and [T5] As shown in Section IIII1 all 
FG-LDPC codes found in [ijj can be constructed using 
our method. The subdesign deletion techniques further 
give infinitely many new codes by fine-tuning their pa- 
rameters and error correction performance. In this sense, 
our method gives many kinds of new and known excellent 
EAQECCs in a single framework. 



V. CONCLUSION 

We have developed a general framework for construct- 
ing entanglement-assisted quantum LDPC codes using 
combinatorial design theory. Our constructions gener- 
ate infinitely many new codes with various desirable 
properties such as high error correction performance, 
high rates, and requiring only one initial entanglement 
bit. Our methods are flexible and allow us to de- 
sign EAQECCs with desirable properties while requir- 



ing prescribed amounts of entanglement. All quantum 
codes constructed in this paper can be efficiently decoded 
through the sum-product algorithm. 

We have introduced many new families of 
entanglement-assisted quantum LDPC codes based 
on combinatorial designs as well as determined all 
fundamental parameters of the well-known families of 
LDPC codes based on finite geometries for most cases. 
Because the entanglement-assisted stabilizer formalism 
bridges classical and quantum codes in a direct manner, 
these results on entanglement-assisted quantum LDPC 
codes are useful both in quantum and classical coding 
theories. 

Our framework encompasses many previously pro- 
posed excellent quantum LDPC codes as well. In fact, 
our method can also be applied to quantum LDPC codes 
under the standard stabilizer formalism by employing the 
ideas found in [HIH]. 

We have focused on the fundamental classes of com- 
binatorial designs. However, other classes of incidence 
structures may provide interesting results as well. For ex- 
ample, the entanglement-assisted quantum LDPC codes 
presented in [l6j can be seen as incidence structures gen- 
erated from the so-called difference matrices and their 
generalizations (see 47] for the definition and basic facts 
about difference matrices) . More general families of com- 
binatorial designs can have nested structures or similar 
strong orthogonal relations between two incidence matri- 
ces. This kind of structure can give asymmetric quantum 
codes (see [57], HH). Structures in finite geometry we did 
not employ may also give interesting quantum LDPC 
codes as well as classical LDPC codes (see, for exam- 
ple, [59|, |60(). Because LDPC codes and sparse incidence 
structures are equivalent, we expect that our methods 
may be further generalized to encompass a wider range 
of both new and known quantum LDPC codes in future 
work. 
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Appendix A: Existence of 2-designs 

Here we discuss the existence of 2-designs to be applied 
to our constructions given in Subsection lllBI The follow- 
ing is the well-known asymptotic existence theorem. 
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Theorem 47 (Wilson [61-63]) The necessary condi- 
tions for the existence of a 2-(v, fi, A) design, \(v — 1) = 
(mod (J, — 1) and \v(v — 1) = (mod — 1)), are also 
sufficient if v > v^x, where v^x is a constant depending 
only on \x and A. 

For fi £ {3,4,5}, necessary and sufficient conditions 
for the existence of an 5(2, /a, v) are known. 

Theorem 48 (Kirkman [HI]) There exists an STS(v) 
if and only if v = 1,3 (mod 6). 

Theorem 49 (Hanani [HI]) There exists an S(2,A,v) 
if and only if v = 1,4 (mod 12). 

Theorem 50 (Hanani [66}) There exists an S(2,5,v) 
if and only if v = 1,5 (mod 20). 

For /i > 6, the necessary and sufficient conditions on 
v for the existence of an S(2, fi, v) are not known in gen- 
eral, although for small values of /j, substantial results 
are known. For a comprehensive table of known Steiner 
2-designs, see [47j . 



Theorems |47l |48l |49l and [50] were proved by con- 
structive methods. Hence, these existence results al- 
low us to construct infinitely many explicit examples of 
entanglement-assisted quantum LDPC codes. It is worth 
mentioning that many of the known proofs of these the- 
orems employ the same construction technique we used 
in Theorem |2"T1 In fact, most S(2,fi, v)s in the original 
proofs of these existence theorems have either Steiner 
spreads or nontrivial subdesigns. 

Numerous other constructions for 2-designs also give 
explicit examples of S(2,/u,,v)s for a wide range of pa- 
rameters. A detailed treatment of STS(v)s is available in 
[67| . Various constructions for S(2, (j,, v)s for many values 
of /j, are also given in [Hj] . 

Appendix B: Parameters of quantum and classical 
FG-LDPC codes with girth six 

Here we give tables of parameters of LDPC codes with 
girth six based on finite geometries. Table IXlVl gives pa- 
rameters of entanglement-assisted quantum LDPC codes 
obtained from PG\(m,q), AGx(m,q), and EGi(m,q). 
Parameters of the corresponding classical FG-LDPC 
codes are listed in Table IXVI 
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TABLE XIV. Parameters of entanglement-assisted quantum LDPC codes from finite geometriesQ 
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a All codes are [[n,fc,d;c]] EAQECCs obtained from PGi(m,q), AGi(m, q), or EGi(m,q). We omit EAQECCs which are created by 
subdesign deletion techniques or do not have dimension greater than one. <fi(m, 2*) is given by Theorem [22] in Subsection IIIIAI 
g(m, 2 t ) is defined as g(m, 2*) = ip(m, 2*) — tp{m — 1, 2*). 

b Type refers to the traditional classification of FG-LDPC codes: Type I uses a line-by-point incidence matrix, while Type II uses the 
transposed (i.e., point-by-line) incidence matrix. 

c The codes obtained from either orientation of the incidence matrix are identical [fSj . 



TABLE XV. Parameters of classical FG-LDPC codesQ 
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a We omit the cases when codes are created by subdesign deletion techniques or do not have enough dimension. 
b The codes obtained from either orientation of the incidence matrix are identical 
c If Hamada's conjecture on the rank of EGi(m, q) [35| is true, the equation holds. 
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